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Dedicated to my Fiance Muhammad Atiq Jamil 

Abstract. The theory of standard bases in polynomial rings with coefRcients 
in a ring R with respect to local orderings is developed, i? is a commutative 
Noetherian ring with 1 and we assume that linear equations are solvable in R. 

1. Introduction 

The aim of this paper is to develop the theory of standard bases especially for 
non-global orderings for polynomial rings with coefficients in a ring. We generalize 
the concept of Adams and Loustaunau ([!]) and Greuel and Pfister ([2]). In the 
book of Adams and Loustaunau the concept of Grobner bases over polynomial rings 
with coefficients in a ring is developed, i.e, they consider standard bases with respect 
to global orderings. In the book of Greuel and Pfister the concept of standard bases 
over polynomial rings with coefficients in a field is developed, i.e, they consider also 
non-global orderings. We will generalize both concepts to a uniform theory. 
Note that the theory of standard bases for ideals developed in this paper can also 
be also developed for modules without changing the proofs. 

First of all we will prove that in the general case the computation of a standard 
basis with respect to a non-global ordering can be reduced using homogenization to 
the computation of a Grobner basis with respect to a suitable global ordering. This 
is also here a very expensive way to compute a standard basis. Therefore later a 
more efficient algorithm similar to [2] is presented. 

Standard basis computations over the rings Z and Z/<m>can be performed using 
the computer algebra system SINGULAR (cf. [4]). 

Standard bases are useful in computing elimination of variables, intersection of 
ideals, quotient of ideals, kernel of the ring map. 

This can be done using the method described in [2] , all the results remain the same 
in our case. 

2. Basic Definitions 

Let i? be a Noetherian commutative ring with 1 and R[xi, . . . the polynomial 
ring in n variables with coefficients in R. Assume that linear equations are solvable 
in R. 
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Definition 2.1. Linear equations are solvable in R if the following conditions are 
satisfied for any a, Oi, . . . , G i?: 

(1) there is an algorithm to compute generators for the R-module 

syzR{ai, ...,am)^ {{h, . . . ,bm) e R"" \ aibi + . . . + a^6^ = 0}, 

(2) there is an algorithm to determine if a & (oi, . . . , am), 

(3) there is an algorithm to compute 6i, . . . , 6^ e such that a — b^ai + . . . + 6^a^ 
if ae (ai, . . . ,a^). 

Wc will use the notations from [1] and [2] and repeat them here for the convenience 

of the reader. 

Definition 2.2. A monomial ordering > is a total ordering on the set of monomials 
MoUn = {x" = Xi^ • . . . ■ x'^"\a = («!, . . . , «„) G N"} in n variables satisfying 

x°' > x^ =^ x'^x'^ > x^'x^ 

for all a, (3, J G N". We also say, > is a monomial ordering on R[xi, . . . , if > 

is a monomial ordering on MoUn- 

Example 2.3. The local lexicographical ordering >is on Mon„ is defined as follows, 
x" > x^ <^==^ 3 1 < i < n, tti = . . . , Q;j_i = ai < Pi. 

Example 2.4. Let M be an invertible {n x n)-matrix with real coefficients and 
Ml, . . . , Mn the rows of M . The matrix M defines an ordering > on Mon„ as 
follows: 

> x^ <^ 31 <i<n, Mia = Mi/5, . . . , M^-io; = Mj_i/3, M^a > Mi(3. 
Every ordering can be defined by a matrix {cf. [2]). 

Definition 2.5. Let > be a fixed monomial ordering. Writing f G R[xi, . . . ,Xn\, 
f ^ 0, in a unique way as a sum of non-zero terms 

f = a„,x"i + a«2^"^ + ■ ■ ■ + aasX'^' , x'^' > x'^^ > . . . > , 

and ttai , ciaa 5 • • • ) (^os ^ ^all: 
LM{f) := x°'^ , the leading monomial of f , 
LE{f ) := ai, the leading exponent of f , 
LT{f) :— ttaiX""^, the leading term of f, 
LC{f) :— ttai, the leading coef ficient of f, 
ecartif) := deg{f) - deg{LM{f)) . 

We define the leading monomial and the leading term of to be 0, and to be 
smaller than any monomial. 

Definition 2.6. Let > be a monomial ordering on Mon„, > is a called global {resp. 
local) ordering if x"" > 1 {resp. < 1) for all (0, . . . , 0). 
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Definition 2.7. For any monomial ordering > on Mon„, 

S> := {u e R[xi, Xn]\{0}\LTiu) = 1} 
is a multiplicatively closed set , 

R[xi, . . .,Xn]> := Sz}R[xi, . . . = {-\f^9 ^ R[xi, ■■■,Xn] andg E S>} 

is the localisation of R[xi, • • • , x^] with respect to Sy and we call R[xi, . . . , 3;„]> the 
ring associated to R[xi, . . . , Xn] and >. 

Definition 2.8. Let > be any monomial ordering. 

For f G R[xi, . . . ,Xn]>, choose u G R[xi, . . . , x„] such that LT{u) = 1 and uf G 
R[xi, . . . ,Xn]- Then 

LM{f) ■.= LM{uf), 
LC{f) ■.= LC{uf), 
LT{f) ■.= LT{uf), 
LE{f) ■.= LE{uf). 

Definition 2.9. Let > be any monomial ordering then, for each G C R[xi, . . . , Xn]> 

LiG) = {{LTig)\geG})n[.,,...,.^] 
is called the leading ideal of G. 

Definition 2.10. Let I < R[xi, . . . , Xn\>- 

(1) . A finite set G C -R[xi, . . . , x„]> is called a standard basis of I with respect to > 
ifGcI, andL{I) =L{G). 

(2) . G is called a strong standard basi^, of I with respect to >, if G <Z I and for 
any f G /\{0} there exists i G {1, . . . ,t} such that LT{gi) divides LT{f). 

(3) . If > is global, a standard basis is also called a Grobner basis. 

3. Computing Standard Bases By Using Homogenization 

Theorem 3.1. Let /i, . . . , G R[xi, . . . , Xn] generating the ideal I < R[xi, . . . , x„]>, 
where > is a monomial ordering given by a matrix M. Let Fi := fl^ G -R[t, xi, . . . , 

be the homogenization of fi and >h be the monomial ordering given by the matrix 
f\ 1 ■■■ 1 
' 

; M 
Vo 

Let {Gi, . . . , Gs} be a Grobner basis, respectively strong Grobner basis of 



""^Strong standard bases do not exist in general. They exist always if i? is a principal ideal 
domain, (cf. theorem 6.4). A strong standard basis is a standard basis. 
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J= {Fi, . . . , Fm) with respect to >h- If we denote Qi = Gj|t=i then, {gi, . . . , gg} is a 
standard basis, respectively strong standard basis of the ideal I with respect to >. 

Proof. Assume that {Gi, . . . , Gg} is a Grobner basis with respect to >h. 

Let / G lnR[xi, . . . , Xn]- Then there exists it e S'> and r}i e R[xi, . . . , Xn] such that 

m 

u - f = ^r]i- fi. 

1=1 

Then there exists p, pi '^l^^p^ Pi>. such that 

m m 

1=1 i=l 

As {Gi, . . . , Gs} is a Grobner basis J so there exist e R[xi, . . . , x„] such that 

s 

LT{tP-u^-f^^) = Y,i,-LT{G,) 

i=l 

putting t — 1 we obtain the result. 

Now assume that {Gi, . . . ,Gs} C J is a strong Grobner basis with respect to > and 
let / G / n R[xi, .... Xn]. We want to show there exists i such that LT{gi)\LT{f) 
and that gi,...,gs& IR[xi, . . . , a;„]>. 
As {Gi,...,Gg}CJ 

m 

with e R[t,Xi, ...,Xn]. 

Put i = 1, we get 

m 

this implies gi,...,gs& IR[xi, . . . , a;„]>. 
Now for / there exists w e S'> such that 

m 

w-f = J2^rfj 

for suitable rjj G it![xi, . . . , Xn]- Then there exists p, pi G Z, p, pj > such that 

m 

i=i 

As {Gi, . . . , G^} is a strong Grobner basis of J there exists i such that 
LT{Gi)\LT{tP ■ ■ f^). This implies LT{gi)\LT{f). □ 
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4. Normal Form 

The concept of a normal form with respect to a given system of polynomials is 
the basis of the theory of standard bases. Normal forms for non-global orderings 
are different and more complicated than normal forms for global orderings. This is 
already the case for polynomial rings over a field. 

Definition 4.1. Let Q denote the set of all finite lists G C R[xi, . . . , Xn]>, 
NF : R[xi, Xn]> R[xi, . . . , x„]>, (/, G) ^ NF{f\G), 

is called a normal form on R[xi, . . . ,Xn]> if for all G d Q and for all 

f e R[xi, . . .,Xn]>, 

(1) . NF{Q\G) = 0, 

(2) . NF{f\G) ^ ^ LT{NF{f\G)) i L{G). 

(3) . If G = {gi, . . . , Qs}, then there exists m e S'> such that r := uf — NF{f\G) 
has a standard representation with respect to G, that is, 

s 

for suitable G R[xi, . . . , Xn] and LM{r) = maxl^^{LM{^i)LM{gi)} . 

To prove the existence of a normal form we give an algorithm to compute it. 

Definition 4.2. Let h G R[xi, . . . ,Xj^, T C . . . , be finite. If h = 0, let 

S{T,h) = 0. Ifh^O, let 

S{T, h) := {EgeT cg ■ • g\ LT{Y.g^^ Cg ■ • g) = LT{h), Cg e R and • LM{cg ■ 
g)^LM{h) ifcg^Q}. 

Renicirk 4.3. The set S{T, h) can be infinite. Algorithm 1 requires to choose an 
element of S{T, h) which is of minimal ecart. This is achieved by computing a gen- 
erating system of S(T, h), which is a kind of a syzygy module, so by the assumption 

on R this can be done. 

Algorithm 1. NF{f\G) 

Let > be any monomial ordering. 

Input: f G R[xi, ...,Xn\, G = {gi, ...,gs} C R[xi, ...,Xn] with gi V i = 
l,...,s. 

Output: h G R[xi, . . . , Xn] a normal form of f with respect to G. 

• h := /; 

• T := G; 

• while{S{T, h) ^ 0) 

choose k G S(T, h) such that ecart{k) is minimal; 
if{ecart{k) > ecart{h)) 
T:=TU{h}; 



6 



AFSHAN SADIQ* 



h := h — k; 

• return h; 

Proposition 4.4. The algorithm terminates and defines a normal form. 

Proof. Termination is proved by using liomogenization with respect to t: 
We start witli h := and T'^ := {g^'lg G G}. 
The while loop looks as follows 

• whi[e(S(T'',t'^h) ^ for some a) 

choose k e S'(T'*, t^h) such that o; > is minimal; 
if(a > 0) 

rph ._ rph y 

h := r/i - k; 

h := {hU=,r, 

Since R[xi. . . . ,Xn]y is Noetherian, there exists some positive integer such that 
L{Tl^) becomes stable for v > N, where denotes the set T'^ after the f-th turn 
of the while loop. The next h, therefore, satisfies LT{h) e L{T^) — L{T'^), whence, 

= LT{h) for some X^ggjife CgX"^ g and a — 0. That is, itself becomes stable for 
V > N and the algorithm continues with fixed T^. Then it terminates, since >h is 
a well ordering on R[t, xi, . . . , Xn]- For the correctness consider the v-th while loop 
of Algorithm 1. There we create (with Hq := /) 

hy = hy-i- y^^CgX^^g 

for some J^geT^g^^^Q ^^^h that LTi^g^r^x"'^ g) = LT{hy-i) and x'^<'LM{cgg) = 
LM{hv-i) if 7^ 0, from the construction of we have 

s v-2 

Cgx'^^g = J2 CiX'^'di + djX^'^o 
geT 1=1 

where 

_ j ai 

which imphes 

s 

hy = - (Y CiX'^'9i + Y djx'^'hj). 

1=1 j=0 

Especially for > 2 we have LM{f) > LM{hy_i) = LM{x'^^Cigi) if Q 7^ and 
LM{f) > LM{h,_i) 



j=0 

9^9i 
g^hj 

9 = 9i 
g^hj 



STANDARD BASES OVER RINGS 



= LM^xf^idjhj) if dj ^ 0. This implies especially x''^^ < 1. 

Suppose by induction, that in the first v — 1 steps (i; > 1) we have constructed 
standard representations 



/ = ^ a^^^Qi + hj, Uj e a'f^ e R[xi, . . . , 



i=l 



with LM{ujf - hj) = maxl^i{LM{a\^^)LM{gi)} ior < j < v - 1, starting with 
uo = 1, af ^ = 0. 

Consider this standard representation ior j — v — 1. 

We replace h^-i by hy + (X]i=i CiX°''9i + Z]J=o djX^^hj), hence we obtain 

s s v—2 

where each hj has a standard representation as above 

s s v—2 s 



=1 i=l j=Q i=l 

v—2 s s v—2 s 

j=0 1=1 i=l j=0 1=1 

s v-2 



iv-1 - Y djX^"^j)f = iY ^i" + Y - Y Yl + K- 



= Y(^t'^ + - E dj^^'<''f^)9i + hv. 

i=l j=Q 

Let Uy := {u^^i - Yl'jZl djxl^^uj) and of-' := af~^^ + Cix"'^ - Yl'jZldjX^'a'f\ We 
have to prove that G and 

s 

Uvf = E ^i^''^' + 
i=l 

is a standard representation, i.e, 

LM(/) = LM{uJ - K) = max|^i{LM(af ^)LM(^i)}. Since a;^ < 1 in case dj ^ 
it follows 

Since LM{f) > LM{x'^'^Ckgk) if 7^ and LM{f) > LM{afgk), x^^ < 1, it follows 
LM{a^^^gk) < LM{f). If LM {a^"'^^ gi) = LM{f) then with the same argument we 
obtain LMia^^^gi) = LM(af "^^^i) = LM(f). □ 

Example 4.5. We consider R = Z and use the local lexicographical ordering Is with 
X > y in Z[a;, y]. Let f — xy^ — 12x'^ then ecart{f) — and let G — {fi, f2, fs} 

where 

fi = -3x + xy, f2 = y^ - 2x^y, /s = <ox^ - x^y^ 
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then ecart{fi) = l,ecart{f2) = l,ecart(/3) = 3. 
In step 1: 

ho = xy^ - Vlx^, T -.^G and 

S(T, ho) = {ky^h + (3A; + V2, k e Z}. 

All elements in S{T, ho) have ecart 1 and we choose xy'^ — xy^ + 4x^y^ e <S'(T, ho). 
Since ecart{ho) < ecart{xy^ — xy^ + Ax^y^) we have to enlarge T: 

T = TU {/4 := ho} and 

hi = xy'^ ~ 12x^ — [xy'^ — xy^ + 4,x^y^) = xy^ — 12x'^ — 4x^y^. 
In step 2: 

hi = xy^ - 12x^ - Ax^y^ with ecart{hi) = 0, 

S{T, hi) = {ky^'h + lxy''f2 + (1 + 3A; - l)yU, k, I G Z}. 

We choose xy^ — Vlx'^y G 5'(T, hi) with minimal ecart and obtain 
h2 = xy^ - Vlx^ - Ax^y^ - {xy^ - 12x^y) = -12x2 + 12x'^y - Ax^y^. 

In step 3: 

h2 = —12x2 + 12x'^y — 4x^y^ with ecart{h2) = 4, 

S{T, h2) = {{2k + 4)xfi + kf2, k G Z}. 

We choose —12x2 ^ 4x'^y G S{T, /i2) with minimal ecart 1 and obtain 
/13 = -12x2 + 12x'^y - Ax^y^ - (-12x2 + Ax'^y) = 8x2y - Ax^y^. 
In step 4- 

hs = Sx'^y - Ax^y^ and LT{h^) i L{T), thus NF{f\G) = Sx'^y - Ax^y^. 

Remark 4.6. Assume R has the following property: c = aiXi + . . . + asXs is solvable 
in R if and only if there exists j and x G i? such that c — UjX. 
Then normal form algorithm is similar to the corresponding normal form algorithm 
for a polynomial ring over a fi,eld, i.e, S{T,h) can be replaced by S{T,h) = {(7 G 
T\ LT{g)\LT{h)}. In this case each standard basis is a strong standard basis. If 
R is a discrete valuation ring or R — Z/(p"^), p a prime number, then R has the 
property above. 

5. Computing Standard Bases 

Theorem 5.1. Let I < R[xi, . . . , x„]> and let G = {gi, . . . , gt} be a set of non-zero 
polynomials in I. Then the following are equivalent. 

(1) . L{G) = L{I). 

(2) . For any polynomial f G i?[xi, . . . , x„]>; we have f E I if and only if 

NF{f\G) = 0. 

(3) . For all f E I , uf = Xlj=i hi9i f^f some polynomials u, hi, . . . , ht G R[xi, . . . , x„] 

such that LT{u) = 1 and 

LM{f) = maxU{LM{K)LM{g,)}. 
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Proof. (1) =^ (2). We know that if NF{f\G) = 0, then f E I. Conversely assume 
that f e I. Let r = NF{f\G) and assume r 0. Since G G I we have r e I. This 
imphes LT(r) e L(I) — L{G). This is a contradiction to Definition 9. 

(2) =^ (3). This is obvious from the Definition of normal form. 

(3) =^ (1). For / G / we need to show that LT(f) G L{G). We have that 
uf = Y^\=ihigi such that LM{f) = maxl^^{LM(hi)LM(gi)}. It is easily seen 
that LT{f) — ^ LT{hi)LT{gi) where the sum is over all i satisfying LM{f) = 
LM{hi)LM{gi). Thus LT{f) e L{G), as desired. □ 

Corollary 5.2. If G is a standard basis of I then I is generated by G. 

Proof. Clearly {gi, . . . ,gt) C /, since each gi is in /. For the other inclusion, let 
f e I. By Theorem 5.1, NF{f\G) = and hence uf G {gi, . . . , gt) R[xi,...,xr.\ for a 
suitable u E □ 

Remark 5.3. Let K be a field and R — K[yi, . . . , y^] be the polynomial ring with 
variables yi, . . . , ym- Let > be a product ordering on K[yi, . . . , ym, Xi, . . . , Xn] such 
that the dominate yi, . . . .y.^ and the restriction of > to R is global. Let 

I < K[yi, . . . , y^n, xi, . . . , Xn]> an ideal and G = {/i, . . . , fs} a standard basis of I. 
Let >i be the ordering on R[xi, . . . ,Xn\ considered as polynomial ring with coeffi- 
cients in R induced by >. Then G is a standard basis of I with respect to >i. 

Proof Let f e I then NF>{f\G) = {NF> the normal form with respect to >). An- 
alyzing the algorithm for iVF> wc obtain Oi, . . . , G /? such that '}2i=i ^i-^^>i(/j) — 
LT^^{f). Here LT^^ is the leading term with respect to the ordering >i. □ 

We use the following Definition from [1]. 

Definition 5.4. Given monomials x"^, . . . , x"* and non-zero elements Ci, . . . , Cg in 

R set L ). Then for a given monomial x'^, we call a syzygy 

h = {hi,...,hs) G syz{L) C (i?[a;i, . . . , a;„])* homogeneous of degree x" provided 
that each hi is a term and x°'^LM{h,i) = for all i such that hi ^ 0. 

Theorem 5.5. Let G = {gi, . . . , gt} be a set of non-zero polynomials in R[xi, . . . , Xn]. 
Let B be a homogeneous generating set for syz{LT{gi), . . . , LT{gt)). Then G is a 
standard basis for the ideal {gi, . . . , gt)R[xi, . . . , x„]> if and only if for all 
{hi, . . . ,ht) G we have 

t 

NF{J2h^9^\G) = 0. 

i=l 

Proof. If G is a standard basis, then by Theorem 5.1, 

t 

NF{J2higi\G)^0. 

i=l 
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Conversely, let e {gi, . . . , gt)R[ ]>, then there exists u E and Ui E 

R[xi, . . . , Xn] such that 

t 

ug = ^Uigi. 

i=l 

Choose a representation as in above equation with = max\^'^[LM{ui)LM[gi)) 
minimal. Since by Theorem 5.1, we need to show that LM{g) = x". We assume 
LM[g) < a;" and show that we can obtain an equation for g with a smaller value 
for x". Let S ^{ie{l,.. . ,t}\LM{ui)LM{gi) = x"}. Then 

Y,LT{ui)LT{gi) = 0. 
ies 

Let h — JZies LT{ui)ei (where ei = (1, . . . , 0), . . . , = (0, . . . , 1) is a generating 
set for R[xx, . . . , Then h e syz{LT{gi), . . . , LT{gt)) and h is homogeneous of 

degree x". Now let B — {hi, . . . , h;}, with hj — {hij, . . . , htj) then h = Yl]=i ^j^j- 
Since h is a homogeneous syzygy, we may assume that the aj are terms such that 
LM{aj)LM{hij)LM{gi) = x" for all i,j such that ajhij ^ 0. By hypothesis, for 
each j, NF(^-^^ hijgi\G) — 0. Thus by Theorem 5.1, for each j = 1, . . . , / there 
exist Wj e S> and Vij e R[xi, . . . , x„] such that 

t 

Wj^hijgi = ^Vijgi, 
ieS i=i 

and 

t 

max\^^LM{vijgi) = LM{^hijgi) < maxl^^LM{hij)LM{gi). 

1=1 

The latter strict inequality is because hijLT{gi) = 0. 

We may assume w — Wj for all j. 

Thus, 

t 

wug = wuigi 

i=l 

= wLT{ui)gi + ^{wui - wLT{ui))gi + wUigi 
ies ies i^s 

I 

— wttjhijgi + terms lower than x^ 
I t 

— Qji'ijS'i + terms lower than x'^. 
j=i i=i 
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Wehaye maXijLM{aj)LM{vij)LM{gi) < maXijLM{aj)LM{hij)LM{gi) = x"". We 
have a representation of as a linear combination of the gi such that the maximum 
of the leading monomials of any summand is less than x". Thus the theorem is 
proved. □ 

As a consequence of Theorem 5.5, we obtain that the following algorithm computes 
a standard basis for a given ideal / in R[xi, . . . , x„]>. 

Algorithm 2. StandardhasisiG) 

Input: F = {/i, . . . , fr} C R[xi, . . . , x„] with fi {1 < i < r), 
Output: G a Standard basis for (/i, . . . , fr)R[xi, . . . , x„]>. 

• G:= F; 

• P:= a finite homogeneous generating set (considered as ordered set) for 

syz{{LT{fi)}i<i<r); 

• while{P ^ 0) 

LetG = {f,,...,fk} 
choose (si, . . . , Sfc ) e P 
P:= P\{(si,...,Sfc)}; 
h:^ NF{Y.UsJ,\G); 
tf{h + 0); 
/fc+i := 

G ■— {/i, . . • , /fc+i}; 

H:— a finite homogeneous generating set for syz{{LT{fi)}i<i<k+i); 
P:^{Px {0}) U {/i = {hi, hk+i) e H\hk+i ^ 0}; 

• return G; 

Using ideas of M.MoUer ([3]) Adams and Loustaunau propose a more efficient 
algorithm for computing Grobner bases. This applies also in our situation with the 
same proof. 

We use the following Definition and Theorem from [1]. 

Definition 5.6. Let x"'^ , . . . , x"'" be a set of monomials. For any subset J C 
{l,...,s}, set x'^-' = lcm{x°''j\j G J). We say that J is saturated with respect to 
x"^ , • • • , x'^" provided that for all j G {1, . . . , s} if x'"-' divides x^'-' , then j G J. For 
any subset J G {!,...,«} we call the saturation of J the set Ji consisting of all 
j G {!,...,«} such that x""^ divides x^-' . (Note that x'^'' — x'^'-'i .) 

Theorem 5.7. (c/. [1], page 214) Given monomials . . . and non-zero ele- 
ments Ci, . . . , Cg in R. For each set J C {1, . . . , s}, which is saturated with respect 
to . . . , x°'% let Bj = • • • , b^^^j} be a set of generators of the R-m,odule of 

syzygies syzji(cj\j G J). (Note that each of the vectors b^j is in the R-module 
where \ J\ denotes the cardinality of J). For each such b^^j, denote its jth coordinate, 
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forj e J, by hjy^j and x^-' = lcm{x°'^\j e J). Set 

x'yj 

( Note that each of the vectors s^^j is in R[xi, . . . ,XnY )• Then the vectors s„_j, for 
J running over all such saturated subsets of {1, ... , s}, and, 1 < i' < I'j, forms a 
homogeneous generating set for the syzygy module syz{cix°''^, . . . , CgX"'"). 

Example 5.8. We consider R = Z and let c\x"'^ = Sxy"^, C2X°'^ = 7xyz, c^x"^ — 
2y^z'^. The saturated subsets of{l, 2, 3} are {1}, {2}, {3}, {1, 2} and {1, 2, 3}. Since 
Z is an integral domain, the singletons {1}, {2}, {3} do not give rise to any non-zero 

syzygy. 

For J = (1, 2} we need to solve in R = Z the equation 3bi + 762 = 0. The module of 
all solutions is generated by (7, —3). Since x^-' — xy'^z, the corresponding syzygy is 

s.,^r^ e,+3^ e,^{-lzM- 

Now for J — {1,2,3} we need to solve 2>bi + 7^2 + 263 = 0. The module of all 
solutions is generated by (—4,2,-1) and (—7,3,0). Then with x''-' — xy^z^ we 
obtain the syzygies are 

^v,j 4 _^^2 6i~l~2 62 y2^^i ^3 ( .^'iyz.i x). 

s.,j = -7^ e, + 3^ 62 = (-7;^^3y;^,0). 
So we obtain that 

syz{3xy\ 7xyz, 2yh') = {{-7z, 3y, 0), {-Az^, 2yz, -x), {-7z\ 3yz, 0)). 
The theorem is the basis of the following modified standard basis algorithm. 
Algorithm 3. Standardbasis(G) 

Input: F = {fi,...,fs}^R[ ] with fiy^O (1 <i <s), 

Output: G a standard basis for (/i, . . . , fs)R[xi.i . . . , a;„]>. 

• G := F; 

• a := 1; 

• m := s; 

• while(a < m) 

Compute -S" = { subsets of {1, ... , a}, saturated with respect to 

LM(/i), . . . , LM{fs), which contain a }; 
for{J e S) 

x^ ■.= lcm{LM{fj)\jeJ); 
Compute a generating set {hi j,i — 1, . . . , /ij} 
for {LC{fj)\j eJ,jyl a)R : (LC(/,))^ 
for(i := 

Compute bj & R, j E J, j ^ a 

such that EjeJj^abjLCifj) + bi,jLC{f,) = 



STANDARD BASES OVER RINGS 13 

ij(r ^ o; 

/m+l '■— f", 

m :— m-\- 1; 

(7 := (7 + 1; 
• return G; 

Example 5.9. We consider i? = Q[a;, y] with the local lexicographical ordering Is in 
R[z\. Let G = {/i, and I = {G) where 

fi=y — x^z^ and /2 = x^y — z 

Then ecart{fi) = 2 and ecart{f2) — 1. 

In step 1 (a = 1) : 

S = {{1}} = saturated subsets of {1} containing 1. 

Since R is a domain and fi , /2 are irreducible we have no non trivial syzygy. 

In step 2 {a = 2) : 

S = {{1.2}} = saturated subsets of {1,2} containing 2. 
J={1,2}. 
x'^ = lcm{l, 1) = 1 

A generating set for (y) : {x^y) is {1}. 

The solution of ybi + x'^y = is bi = —x^. 
h = —x'^{y — x^z'^) + l{x'^y — z) = —z + x^z"^ 
which is reduced with respect to G = {/i, /2}- 
/3 h, 

and G = {/i,/2,/3}- 

In step 3 ia = 3) : 

S = {{1, 2, 3}} = saturated subsets of {1, 2, 3} containing 3. 
J ={1,2, 3}. 
x'^ — lcm{l, 1, z) = z. 

A generating set for {y^x'^y) : ( — 1) is {y}. 

The solution of ybi + x^yb2 + x^y = is bi = l,b2 = 0. 
h — z{y — x^z^) + y{—z + x^z"^) = x^yz^ — x^z^. 
NF{h\G) = 0. 

So G — {/i, /2, fs} is a standard basis. 

Remark 5.10. Let R be a local ring of the type K[yi, . . . ,ym]{yi,...,y^)/I , where 
I <K\ Xi, ■ ■ ■ , Xn]{xi,...,xn)- Then we can compute Hilbert- Samuel function of ideals 
in R[xi, . . . ,Xn]> by using the same method described in [2]. Carlo Traversa [5], 
used the Hilbert- Samuel function to speed up the Buchberger Algorithm. The same 
method can be applied in our case too. 



h := spoly{fi, /s) = -^,{-^y + xy) - —{y"^ - 2x) = 6x - xy' 
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6. Standard Bases Over Principal Ideal Domains 

If i? is a principal ideal domain (for short PID) then there is a standard basis 
algorithm similar to the corresponding algorithm for a polynomial ring over a field 
(cf. [2], page 54) with the following notion of the s-polynomial. 

Definition 6.1. Let f.gE R[xi, . . . , Xn]\{0}. 

lcm{LT{f),LT{g)) ^ lcm{LT{f), LTjg)) 
spoly(f,9) jj^^ / g. 

This is a consequence of [1], Proposition 4.5.3 and Theorem 5.5. 

Example 6.2. We consider R = Z with the local lexicographical ordering Is with 
X > y in 1.[x, y] . 

Let / = (/i, /2) where fi = -3y + xy and /a = _ 2x. 
Initialization: G = {/i, /2}; ^ = {{/i, /2}}- 
In step 1: 

LT{h) = -3y, LT{h) = y\ 

c = Zcm(ci, C2) = /cm(3, 1) = 3, = lcm{LM{fi), LM(/2)) = lcm{y, y^) = y^ 

— {-3y + xy) 

-3?/ y^ 

which is reduced with respect to G. 

SoGU {/3 h} and G = {{/i, /a}, {/2, fs}}. 

In step 2: 

LT{h) = -3y,LT{f,) = 6x, 

c — lcm{ci, C3) = lcm{3, 6) — 6, x'^ — lcm{LM{fi), LM{fz)) — lcm{y, x) — xy 

spolyifi, fs) = 3^(-3|/ + xy) - ^i^x - xy"^) = xy^ - 2x^y = xy ■ f2 

NF{xy^ - 2x^y\G) = 0. 
In step 3: 

LT{h)=y\LT{h) = %x, 

c = lcm{c2, C3) = /cm(l, 6) = 6, = lcm{LM{f2), LM{fs)) = Icm^y"^, x) = xy"^ 

spoly{f2: h) = — - 2x) - — — {6x - xy'^) = xy^ - 12x^ 

ox 

NF{xy^ - 12x'^\G) = 0. Since xy^ - Ux"^ = xy'^f2 - 12x^ + 2x'^y^ = xy^f2 + 2xh. 
So G — {/i, /2, /a} is a standard basis. 

Example 6.3. We consider R = Z with the local degree lexicographical ordering ds 

(c/.[2]) with X > y > z in Z[a;, y, z] . 

Let I = (/i, /2, /a, h) where fi = 15x^ + 28y'^z'^ and f2 = 3x^y + lyz^ /a = Axy^ - 
hxz^^, fi = -28y^ + 35yz^\ 
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Similar to Example 6.2 we have to compute the normal form of the spoly's of all 
pairs ifijj). 

Nf(spoly(fi, f2)\G) = 35yz' - 28yh^ =:f,,G^GU {f,} 

Nf{spoly{f2, /3)|G) = 28y^z' + 15x'z'^ =:f,,G = G\J {f,} 
Nf{spoly{f,j5)\G) = 35x^15 ^. q = GU {M- 
The normal form of all the other spoly's is zero. 

The standard basis is G — {/i, /2, /s, /4, /s, /e^ /?}• G is a standard basis of I in 

Z[a;, y, z] {x,y,z)- 

Note that fi — zf^ = 15x^ — 15x'^yz^^ = 15x^(1 — yz^'^). Therefore 15x'^ G /. 
Similarly one can see that 28y'^ , 28y^ z^ , 35y z^ and 35x2;^^ G /. This implies that 
{15x^, Sx'^y + 7yz^, Axy"^ — bxz^^ , 28y^, 28y'^z^, 35yz^, 35x2;^^} is a standard basis of 
I. 

Theorem 6.4. (c/. [1] Theorem 4.5.9 page 251 ) Let R be a PID, and I be an 

ideal of R[xi, . . . ,a;„]>. Assume that {/i, . . . , /„} is a standard basis for I. Let 
LT{fi) = CiX°'' , for a saturated subset J of {1, . . . , s}, let cj = gcd{cj\j G J) and 
write cj = J^jej'^j'^i {any such representation will do). Also, let x"^-^ = lcm{x°'^\j G 
J). Then the set 

aj—^fj\J is a saturated subset o/ {1, . . . , s}} 

is a strong standard basis for I. In particular, every non-zero ideal in R[xi, . . . , x„]> 
has a strong standard basis. 

7. Standard Bases In The Formal Power Series Rings 

Finally we want to apply our results to compute standard bases in the formal 
power series ring . . . , Xn]] with coefficients in a ring R. Let > be a local degree 

ordering, i.e, > is a local ordering and x°' > x^ implies deg{x°') < deg{xf^). A non- 
zero element / G i?[[xi, . . . , can be written as X^^o (^v^'^'" , & R, ao ^ and 
x"" > x^^+i for all V. As in Definition 2.5, wc define LM{f), LE{f), LT{f), LC{f) 
and tail{f ). As in Definition 2.11, we define a standard basis (respectively a strong 
standard basis) of an ideal / C i?[[xi, . . . , x„]]. 

Proposition 7.1. Let I < R[xi, . . . ,x„] and G is a standard basis (respectively a 
strong standard basis) of I with respect to >, where > is a local degree orderig. Then 
G is a standard basis (respectively a strong standard basis) of IR[[xi, . . . ,x„]]. 

Proof. Let {gi, . . . , gg} be a standard basis of / and g = ^^^^ Oigi G IR[[xi, . . . , x„]], 
g 0. Let c be an integer such that LM{g) ^ (xi, . . . , x^)'^. Choose G 
R[xi, . . . , Xn] such that lii — ai G (xi, . . . , XnY- Let g = J2i=i '^idi- Then g E I 
and g — g E (xi, .... This implies LT{g) = LT{g). If G is a strong standard 

basis for / then there exists i such that LT{gi)\LT{g) = LT{g), i.e, G is also a 
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strong standard basis of IR[[xi, . . . , Similarly it follows that a standard basis 
of / is a standard basis of IR[[xi, . . . , Xn]]- □ 

8. Procedures 

Let / < R[xi, . . . ,Xn] be an ideal and G — {gi, . . . , g^} be a standard basis of /. 
Then we can compute a strong standard basis of / using the SINGULAR-procedures 
below. 

LIB"poly.lib"; 

proc powerSet (int n) / /computes the set of all subsets of {1, ... , n} 
{ 

list L,K,S; 
int i; 
if (n==0) 

{ 

L[1]=L; 
return (L) ; 

} 

if (n==l) 
{ 

L[1]=L; 
L[2]=list(l) ; 
return (L) ; 

} 

S=powerSet (n-1) ; 
int r=size(S) ; 
S[r+l]=list(n) ; 
for(i=2;i<=r;i++) 

{ 

K=S[i] ; 

K[size(K)+l]=n; 
S[size(S)+l]=K; 

} 

return (S) ; 

} 

proc satt (ideal I) / /computes the saturated subsets of {1, ... , size{I)} 

{ //w.r.t {LM{f) I / e /} 

int j; 
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list K,J; 
I=lead(I) ; 

J=powerSet (size(I) ) ; 
f or(j=2; j<=size(J) ; j++) 

{ 

if (saturatCl, J[j] )) 

{ 

K[size(K)+l]=J[j] ; 

} 

} 

return (K) ; 



proc specialGCD2(int a, int n) 
{ 

int x=a mod n; 

if (x==0)return(list(0,l,n)) ; 
list L=specialGCD2(n,x) ; 

return(list(L[2] ,L[l]-(a-x)*L[2] /n,L[3])) ; 



proc specialGCD(list L) //computes the gcd over Z 
{ 

int i; 

for(i=l;i<=size(L) ;i++)L[i]=int(L[i] ) ; 
if (size(L)==l)return(L) ; 

if (size(L)==2)return(specialGCD2(L[l] ,L[2] )) ; 

bigint p=L[size(L)] ; 

L=delete (L , size (L) ) ; 

list T=specialGCD(L) ; 

list S=specialGCD2(T[size(T)] ,p) ; 

f or (i=l ; i<=size (T) -1 ; i++) 

{ 

T[i]=T[i]*S[l] ; 

} 

p=T[size(T)] ; 
T[size(T)]=S [2] ; 
T[size(T)+l]=S[3] ; 
return (T) ; 
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proc coeffJ(list C,list J) // computes the coefficients such that 
{ // 9cd{cj\j e J) = Ejgj ajCj 

if (size(J)==l)return(list(l)) ; 

int n=size(J) ; 

int i; 

list L,M; 

for(i=l;i<=n;i++) 

{ 

L[size(L)+l]=C[J[i]] ; 

} 

L=specialGCD(L) ; 
L=delete(L,size(L)) ; 
return (L) ; 

} 



proc IcmJ (ideal X,list J) 
{ 

poly p=X[J[l]] ; 
int i; 

for(i=2;i<=size(J) ;i++) 
{ 

p=lcmS(p,X[J[i]]); 

} 

return (p) ; 

} 

proc maxZ(intvec a.intvec b) 
{ 

int i,j; 
intvec c; 

f or (i=l ; i<=size(a) ; i++) 
{ 

if (a[i]>=b[i]) 
{ 

c[i]=a[i] ; 

} 

else 
{ 

c[i]=b[i] ; 



/ / computes lcm{fj \ fj & X, j & J) 



1 1 computes the maximum of two 
/ / integer vectors 
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} 

} 

return (c) ; 

} 

proc IcmSCpoly p,poly q) // computes the LCM of two monomials 

intvec a=leadexp(p) ; 
intvec b=leadexp(q) ; 
intvec c=maxZ(a,b); 
int i,j; 
poly s=l; 
f or(i=l ; i<=nvars (basering) ; i++) 

{ 

s=s*var (i) "c [i] ; 

} 

return (s) ; 

} 

proc leadTermC ideal I) // give a list containing two lists, first 

{ // contains leading coefficients and second contains 

int i ; / /leading monomials 

list L; 

ideal J; 

I=lead(l) ; 

f or(i=l ; i<=size(I) ; i++) 
{ 

L[size(L)+l]=leadcoef (I [i] ) ; 
J [size(J)+l] =leadmonom(I [i] ) ; 

} 

returndist (L, J) ) ; 

} 

proc clean (ideal G) // delete from an ideal /j such that 

{ // there exists fj e I with LM{fj) divides LM{fi) 

int i , j ; 

while (i<=size (G) -1) 
{ 

i++; 

while (j < size(G)) 



20 



AFSHAN SADIQ* 



{ 

if ((G[i] !=0)&&(G[j] !=0)) 
{ 

if (lead(G[j])/lead(G[i]) !=0) 
{ 

G[j]=0; 

G=siinplify(G,2) ; 

j--; 

} 

if (G[j] !=0) 
{ 

if (lead(G[i])/lead(G[j]) !=0) 
{ 

G[i]=0; 

G=simplify(G,2) ; 

i— ; 
break; 

} 

} 

} 

} 

} 

return (simplify(G, 2)) ; 

} 

proc StrongSB (ideal I) // computes a strong standard basis 

{ //for an ideal if standard basis is given 

def R=basering; 

list rl=ringlist (R) ; 

rl[l]=0; 

def S=ring(rl) ; 

setring S; 

ideal I=imap(R,I); 

list L=leadTerin(I) ; 

list C=L[1] ; 

ideal X=L [2] ; 

list J=satt(X); 

int i , j ; 

ideal G; 
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list P; 
poly f,q; 

f or(i=l ; i<=size(J) ; i++) 
{ 

P=coeff J(C,J[i]) ; 
q=lcmJ(X,J[i]) ; 
f=0; 

for(j=l; j<=size(J[i]) ; j++) 
{ 

f =f +q/X [J [i] [j] ] *I [J [i] [j] ] *P [j] ; 

} 

G[size(G)+l]=f ; 

} 

setring R; 

ideal G=imap(S,G); 

return (clean (G) ) ; 



Example 8.1. Consider Example 6.3, we now compute strong standard basis for 
I = {15a;2, Sx'^y + 7yz^ Axy^ - Bxz^'^, 28y^ 28y^z^, 35y^^ 35xz^^}. 

> ring R = integer, {x, y, z), ds; 

> ideal I = 15x2, 3x2y + 7yz5, Axy2 - 5x^10, 287/3, 2%y2zh, 35yz5, 35xzl5; 

> strong SB {I); 
[1] = 15x2 

[2] = 3x2y + 7yz5 

[3] = 4xy2 - 5xzl0 

[4] = x2y2 - 7y2z5 - 5x2^10 

[5] = 28y3 

[6] = 35yz5 

[7] = x2yzb + 84y^l0 

[8] = 7y2z5 

[9] = 35xzl5 

[10] = hx2zlh 
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